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Abstract 

We consider the response of a dynamical system driven by external adiabatic 
fluctuations. Based on the 'adiabatic following approximation' we have made a 
systematic separation of time-scales to carry out an expansion in a|/i| _1 , where 
a is the strength of fluctuations and |/x| is the damping rate. We show that 
probability distribution functions obey the differential equations of motion which 
contain third order terms ( beyond the usual Fokker-Planck terms ) leading to 
non-Gaussian noise. The problem of adiabatic fluctuations in velocity space 
which is the counterpart of Brownian motion for fast fluctuations, has been solved 
exactly. The characteristic function and the associated probability distribution 
function are shown to be of stable form. The linear dissipation leads to a steady 
state which is stable and the variances and higher moments are shown to be 
finite. 



PACS NO. : 05.20 Dd 



1 



I. Introduction 



In this paper we have discussed the stochastic dynamics of a system driven by ex- 
ternal, adiabatic fluctuations. The opposite counterpart of these processes correspond 
to stochastic processes with fast fluctuations which are more frequently encountered in 
physical and chemical sciences. The classic and celebrated problem of the latter kind 
is the century-old paradigm of Brownian motion first correctly formulated by Einstein 
[1,2]. In dealing with fast stochastic processes one essentially examines the average 
motion of the system subjected to fast fluctuations ( which may be of external or of 
internal type ) with the following separation of time-scales in mind. If r c is the corre- 
lation time of fluctuations which is the shortest timescale in the dynamics, compared 
to coarse-grained timescale At over which one follows the average evolution, then 



where refers to the inverse of the damping rate ( or inverse of the largest eigenvalue 
of the "unperturbed" system ). Herein we analyze the average dynamics of a general 
multivariate nonlinear system subjected to external, adiabatically slow fluctuations. 
We have derived the equation of motion for evolution of the probability distribution 
function in phase space on a coarse-grained timescale At assuming that At satisfies 
the following inequality 



r c < At < — , 



(i) 



-- < At < r c . 



(II) 
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The slow fluctuations characterized by very long correlation time have received a lot 
of attention of various workers over the years [2,3]. While the treatment of stochastic 
differential equations with fast fluctuation is based on the assumption that there is a 
very short correlation time such that one is allowed to make an appropriate expansion 
in ar c , where a is the strength of fluctuation, simplified assumption for dealing with 
long correlation time is rather relatively scarce. In general, the problem of long cor- 
relation time is handled at the expense of severe restriction on the type of stochastic 
behavior. For instance, several authors [2,3] have tried the linear and nonlinear mod- 
els within the framework of simple Markov processes of the type, dichotomic process, 
two-state Markov process, random telegraphic process, etc. Our strategy here is to 
follow a perturbative approach, pertaining to the separation of the timescale (II) with- 
out keeping any above-mentioned restriction on the type of stochastic behavior. Based 
on the 'adiabatic following approximation' [4] we have recently [5] carried out an ex- 
pansion in to obtain a linear differential equation for the average solution. In 
this paper we extend this analysis to treat nonlinear stochastic differential equations 
for construction of appropriate master equations. The perturbative expansion is es- 
sentially a counterpart of expansion in ar c as dealt in the case of fast fluctuations [2]. 
The difference between the two expansion schemes lies in identification of two distinct 
shortest time-scales in the dynamics of the two cases. In the case of fast fluctuations 
it is t c , whereas the corresponding role is played by in adiabatic fluctuations. 

We have shown that the equation of motion in phase space for probability distribu- 
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tion function contains beyond the ordinary Fokker-Planck terms, third order derivative 
terms. As shown by Pawula [6] for one dimensional case, an equation with third order 
derivative terms is in contradiction to the positivity for transition probability for short 
time. However it is well known that finite derivative terms of order larger than two may 
be quite useful in different occasions [7-8], e.g., in the treatment of optical bistability 
described in terms of quasi-distribution function of Wigner in quantum optics [7] or 
explaining trimolecular reactions using Poisson representation of Fokker-Planck equa- 
tion, also in one dimensional random walk with boundary within a scheme of expansion 
of master equation. Although at this stage of development a clear general probabilistic 
interpretation in terms of any real stochastic process is lacking [7] one can identify 
the noise terms with distinct characteristics for such processes. In a similar spirit we 
are led to the conclusion in the present context that adiabatic fluctuations give rise to 
third order non-Gaussian noise terms beyond the usual Fokker-Planck terms. 

The central result of this paper is the solution of the problem of adiabatic fluc- 
tuations in velocity space, which is the counterpart of Brownian dynamics for rapid 
fluctuations. We have shown that the characteristic function obeys a simple third order 
differential equation. This can be solved exactly to obtain a probability distribution 
of stable form which for small arguments displays a power law behavior. It is also 
important to note that the linear dissipation leads to a stable steady state distribu- 
tion. However, the fluctuation being external the energy supplied by this cannot be 
balanced by dissipation and as such there is no fluctuation-dissipation relation in this 

4 



case. Furthermore, the non-Gaussian statistical characteristics can be obtained from 
the calculation of variances and higher moments which are shown to be finite. We thus 
conclude that although in many cases third order noise makes the probabilistic consid- 
eration truly difficult, the systems driven by adiabatic fluctuations display a distinct 
non-Gaussian stochastic behavior is amenable to understanding in simple probabilistic 
terms. Occasionally wherever possible we allow ourselves a fair comparison with Levy 
processes [9,11] and point out the essential differences. 

The outlay of the paper is as follows : In the next section we review the basic 
aspects of adiabatic fluctuations in linear processes as dealt in our earlier paper [5]. 
The two basic assumptions, e.g., the adiabatic following approximation and the de- 
coupling approximation as well as validity and convergence of perturbative expansion 
were discussed in detail in the earlier paper [5]. To make this paper self-contained and 
readable we review its salient features. In Sec. Ill we extend the treatment to nonlinear 
equations. The equations in phase space have been derived in Sec. IV. The counter- 
part of Brownian motion in velocity space for slow fluctuations have been treated in 
Sec.V. Explicit solution for the probability distribution function and the approach to 
equilibrium have been discussed. The paper is concluded in Sec. VI. 
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II. Linear processes with adiabatic fluctuations 



To begin with we have considered the following linear equation, 

it = {A + aAi(t)}u, (1) 

where u is a vector with n components, A is a constant matrix of dimension n x n 
with negative real eigenvalues and Ai(t) is a random matrix, a is a parameter which 
measures the strength of fluctuation. 

It is convenient to assume that Ai(t) is a stationary process with (Ai(t)) = 0. Eq.(l) 
sets the two time scales of the system, measured by the inverse of the largest eigenvalue 
of the matrix A and the time scale of fluctuations of Ai(t) (more precisely correlation 
time of fluctuation). In the problem of Brownian motion where one deals with very 
fast fluctuations such that correlation time r c is essentially the shortest time scale in 
the dynamics, one thus follows the evolution of the average (u) on a coarse-grained 
timescale. 

Before proceeding further we now make two remarks at this stage : First, since in 
the present context we have considered a stochastic process in which the fluctuations 
are weak and adiabatically slow, Ai(t) is an adiabatic stochastic process. Therefore 
the usual procedure of systematic expansion in ar c which relies on smallness of r c , is 
not valid. We thus take resort to a different approach based on an expansion in a|/x| _1 , 
where refers to the largest eigenvalue of A matrix. Second, we do not make any 
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a priori assumption about the nature of the stochastic process, such as, Gaussian or 
dichotomic etc. The only assumption that have been made about the stochastic process 
is that the inverse of the damping rate is much shorter compared to the correlation 
time of fluctuations Ai(t). 

As a first step we introduce an interaction representation as given by, 

u(t) = exp(A t)v(t) , (2) 
and applying it to Eq.(l) we obtain, 

v = aV(t)v, (3) 

where, 

V(f) = exp(-A„t)Ai(f) exp(A t) . (4) 

On integration Eq.(3) yields, 

v (t) = v (0) + a /* V(t')v(t')dt' . (5) 
Jo 

On iterating Eq.(5) once, we are led to an ensemble average equation of the form, 

(v(t)) = v(0) + a 2 f alt' f dt"(V(t')V(t")v(t")) . (6) 
jo Jo 

The equation is still exact since no second order approximation has been used. 

Now taking the time derivative of v(t) we arrive at the following integrodifferential 
equation in which the initial value v(0) no longer appears, 

j t (v(t)) = a 2 j\vm(t'Ht'))dt'. (7) 



Making use of a change of integration variable t' — t — r and reverting back to the 
original representation we obtain 

4;{u(t)) = A (u) + a 2 /'<Ai(*) exp(A r)A 1 (t - r)u(t - r))dr . (8) 
at Jo 

The adiabatic following assumption (see the discussion at the end of this section), 
that Ax(t) and the components of u(t) vary slowly on the scale of inverse of A , 
can now be utilized. Following Crisp [4] we note that a Taylor series expansion of 
Ai(t — r)u(t — t) in the average (. . .) of the a 2 -term in Eq.(8) allows us to reduce the 
above equation to the following form, 

|(^)) = A ( M ) + a 2 f:^(A 1 (t)I n ^[A 1 (tMt)]) . (9) 
I n can also be written as 

JO ■ 

= E D V -n+1 D Jk , Re »33 < • 

3 % 

Eq.(9) can then be rewritten in the form 

d 00 d n 

-<«(*)> = Ao(«> + o? E(-l) n (Ai(t)DE n+1 D^ — \A x (t)u(t)\) , (10) 
where we use 

I n = n!DE n+1 D" 1 . (11) 
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Here D is a matrix which diagonalises A and 



E 



n+l 



W 1 J 



and fj,jj are the eigenvalues of A . 

Although the Eq.(lO) involves an infinite series it is expected to yield useful re- 
sult in the adiabatic following approximation. Under this approximation the quantity 
[Ai(t)u(t)] varies very little (such that ^(A^) in Eq.(lO) is small) and also since \fj,jj\ 
in E n+ i is large the series in Eq.(lO) ( which is thus an expansion in ) converges 

rapidly. Keeping only the two lowest order terms we arrive at, 

d 



i-Mt)) = A (u) + a 2 (Ai(i)XiAi(i)u(t)) - a 2 (Ai(t)X 2 Ai(t)u(i)> 
at 

-a 2 {A 1 (t)X 2 A 1 (t)u(t)) 



;i2) 



where, 



D E„ . i D 



n+l 



It is evident that the average (it) is related to a more complicated average. Following 
Bourret [10], we now implement the decoupling approximation. This allows us to break 
up the average as a product of averages. Keeping terms only of the order of a 2 we 
obtain 
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d 
di 



(u(t)) = {A + ^[(Ax^XxAxft)) - (A 1 (OX 2 A 1 (i)) 



(A 1 (t)X 2 A 1 (*))A ]} («(*)>. 



(13) 



Thus the average of u(t) obeys a nonstochastic differential equation in which the 
effect of weak adiabatic fluctuations is accounted for by renormalizing A through the 
addition of constant terms of the order of a 2 . 

The implementation of Bourret's decoupling approximation [10] is a major step for 
almost any treatment of multiplicative noise upto date [2,3,12]. This is because of the 
fact that the average of a product of stochastic quantities does not factorize into the 
product of averages, although it has been argued that good approximations can be 
derived by assuming such factorization. In the case of fast fluctuations it has been 
justified if the driving stochastic noise has a short correlation time such that Kubo 
number a 2 r c is very small in the cummulant expansion scheme ( an expansion in ar c 
). The factorization has been shown to be exact in the limit of zero correlation time 
and in some cases of specific noise processes [3,12] and the solution for the average can 
be found exactly. 

In contrast to cummulant expansion (valid in the case of fast fluctuation which relies 
on an expansion in ar c ) the present scheme of adiabatic following approximation results 
in a perturbation series, where the n-th term is of order a-jj^[A(t)u(t)]/ fiff 1 and the 
convergence of the series is assured since the numerator varies little in the scale of 
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Eq.(13) is a result of decoupling approximation employed in this expansion 
scheme. If one neglects the free motion due to Ao term then Eq.(13), which gives the 
lowest order evolution, asserts that 



The contribution of is derived from X x of the first term in Eq.(13), (i.e., due 
to E n+ i matrix). Note that because of full integration over r in moving from Eq.(8) 
to Eq.(9) correlation time r c does not appear in Eq.(13) and the time-scale set by the 
dynamics is l/i]" 1 only. For a fast process on the other hand the counterpart of the last 
relation is [12] 



It is also easy to calculate the relative error made in the decoupling approximation. 
We first note that Eq.(13) is obtained from Eq.(8). To the second order it means 
omitting terms of the order (a At) 3 and higher (where At is the coarse-grained time- 
scale over which (u) evolves). As the lower bound of At is determined by \fi\~ 1 , it 
implies that we neglect terms of the order (a|/i| -1 ) 3 in the evolution equation. The 
relative error made in the decoupling approximation is thus of the order (o^l -1 ) 3 
which is well within the order of lowest order evolution. We thus see that the adiabatic 
expansion is an expansion in and the decoupling approximation in the slow 

fluctuation is valid where a 2 |/z| _1 is very small. Thus u(t) in the average ( in the right 



It 





a 2 r c (u) . 
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hand side of Eq.(12) ) is realized as an average (u(t)) ( which varies in the coarse- 
grained timescale At ) in Eq.(13) pertaining to the separation of the time-scales in the 
inequality (II) in Sec. I. 

Before closing this section a few pertinent points regarding the notion of "adiabatic 
following approximation" and its genesis may be noted. The notion has acquired special 
relevance in the quantum optical context where one is concerned with a two-level atom 
interacting with single mode electromagnetic field. The model is described by the 
standard Bloch equations, where the field strength varies slowly on the time-scale of 
inverse of the damping constant or the frequency detuning between the atom and the 
field. If the field is varying adiabatically enough, then the population inversion of 
the Bloch vector components would follow the field adiabatically in going from ground 
to upper state, i.e., the ground state population is adiabatically inverted. The term 
li adiabatic following" is thus used to describe collectively the associated experimental 
phenomena [19]. 
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III. Probabilistic considerations : Extension to nonlinear 



equations 



We now generalize Eq.(l) to a stochastic nonlinear differential equation written in 
the following form 

Uu = F v {{u v }, t; £(*)) , v = 1,2,..., N . (14) 

The above equation determines a stochastic process with some initial given con- 
dition {u„(0)}. £(t) is the adiabatic stochastic process. It may be pointed out that 
the treatment given in the last section cannot be extended directly to this equation 
to obtain an equation for average (u) since nonlinearity in Eq.(14) results in higher 
moments. However, it is possible to transform the nonlinear problem to a linear one if 
one considers the motion of a representative point-w in n-dimensional space (ui . . . u n ) 
as governed by Eq.(14). The equation of continuity , which expresses the conservation 
of points determines the variation of density in time, 

9j ^f 1 = ~ E ^(K>, t; m)p(u, t) (15) 

or more compactly 

| = -V-F,. (16) 

Eq.(15) is a linear stochastic differential equation and is an ideal candidate for the 
method discussed in the last section for the linear case. We emphasize here that the 
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basis of the present analysis is essentially the two approximations as introduced earlier 
and no further approximation is needed to extend the analysis to nonlinear domain. 
F can now be split as 

F(K},i;e(*)) = Fo(K})+aF 1 (K},t;e(*)) , (17) 

where ¥q{{u u }) is the constant part and Fi({u u }, t; £(£)) is the random part with 
(Fi(£)) = 0; a is the parameter defined earlier which measures the strength of fluctu- 
ation. Eq.(16) therefore takes the following form, 

p(u, t) = (A + aAi)p(«, t) , (18) 

where A = —V ■ F and Ai = —V ■ Fi. The symbol V is used for the operator that 
differentiates everything that comes after it with respect to u. 

With the above identification of A and Ai we are now in a position to apply the 
fundamental Eq.(9) derived in the earlier section, to Eq.(18). We have 



oo (_-\\n An 

-V • F P + a 2 £ ^-tM-V • FJn— (-V • FxP)> 



(19) 



where (p(u,t)) = P{u,t) and also 

POO 

/ dre- rV - Fo r n . (20) 
Jo 



/•oo 



Adiabatic following approximation may now be invoked again in the spirit of earlier 
treatment in Sec. II to obtain 
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F + a 2 (F!l V • Fi) - a 2 (F 1 I 1 V • F^ 



+ a^FxIxV • FxV • F ) P(u,t) , 



(21) 



where we keep terms of the order of a 2 for n = and 1 of the series in Eq.(19). 



Our next task is to simplify further the expressions for the averages in Eq.(21). To 
this end we first note that the operator exp(— rV • F ) provides the solution of the 



equation 



df(u,t) 
dt 



-V-Fo/Kt) , 



(22) 



(/ signifies the unperturbed part of P) which can be found explicitly in terms of 
characteristic curves. The equation 

u = F (u) 

for fixed t determines a mapping from u{r = 0) to u(t), i.e., u — > u T with inverse 
(u T )- T = u. The solution of Eq.(22) is 

d(u"*) 



=/(«"*, 0) 



d(u) 



e- tVF °f(u,0) , 



(23) 



din-*) 



d(u) 



being a Jacobian determinant. The effect of exp(— tV ■ F ) or /(u) is as follows 



exp(-iV.F o )/(u,0) = /(«"*, 0) 



d(u) 



(24) 



The relation (24) may be used to simplify the average in Eq.(21). We thus have 

du~ T 



roc 

(V • FiloV • Fx) = V • / (F!V_ r • F^O) 

Jo 



du 



dr , (25) 
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roc 

(V • FxIxV • Fi) = V • / r(F 1 V_ T • F^O) 

Jo 



du 



du 



dr , (26) 



(V • F1I1V • FiV • F ) = V • / r(F 1 V_ T • Fi( U - T )V_ T • F ( u - T )) 

Jo 



du 



du 



dr . (27) 



The use of Eq.(25)-(27) reduces Eq.(21) to a more tractable form, 



- V-|F -a 2 ^ 00 (F 1 V- r -F 1 (O) 



du T 



du 



dr 



roc 

+ a 2 / r(F 1 V_ T -F 1 ( M - T )) 
Jo 



du 



du 



dr 



a: 



roc 

' r(F 1 V_ r -F 1 ( M - r )V- r -Fo( M - r )) 
Jo 



du 



du 



dr}P(u,t), (28) 



where V_ T denotes the differential with respect to U- T . Eq.(28) is our basic result in 
this section. The equation is second order in a, i.e., of the order of a 2 |/i| -1 , where 
refers to the eigenvalue of A . In our earlier communication [5] we have shown 
the convergence of the series in pertaining to the separation of the time-scales 

implied in II in Sec. I. We also remark that it is possible to extend the treatment 
to higher order, in general. It is also to be noted that the equation involves three 
differentiation of P(u,t) with respect to the components of u and is a third order 
equation. The appearance of third-order noise beyond the usual Fokker-Planck terms 
is a characteristic of the process we consider here. We discuss this aspect in more detail 
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in the following two sections. 



IV. Adiabatic stochasticity in phase space 



We now consider the motion of a particle in one dimension subjected to a force K(x) 
depending on the position x, a frictional force —f3x and a stochastic force ee£(t). Here 
(3 is a measure of damping of the system and a is the strength of adiabatically slow 
fluctuations We thus write 



mi + (3x = K(x) + cx£(t) . 



(29) 



The corresponding problem of fast fluctuation a£(t) was studied by Kramers [13] as 
a model of simple chemical reactions and by Bixon and Zwanzig [14] as a model for 
fluctuating nonlinear systems. 

For simplicity we set m — 1 and x — v. Then the two components of u in this 
example are x and v. Taking Eq.(17) into account we have 

F 0x =v F lx = 

F 0v = -(3v + K(x) F lv = a£(t) 

By considering a small variation of v in time r, one obtain (from the unperturbed 
version of Eq.(29)) the Jacobian determinant for the "unperturbed" mapping u — > u T 



(30) 



du T 




d(x T ,v T ) 


du 




d(x, v) 



1 + (3t + 0(t 2 ) 



(31) 
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and 



ox Ox ox ov v / 

The Eq.(28) now reduces to the following form 



(32) 



d_ 

dt 



P( x , v ,t) = -V • |f - a 2 jH(FiV_ T • Fi(aT T , O) 



+a 2 / r(F 1 V_ r -F 1 (a;- T ,0) 



d(x~ T , f ~ r ) 



/■OO 

-a 2 / r(F 1 V- r -F 1 (x- r ,OV-.-Fo(x- T ,0) 
Jo 



d(x T ,v T ) 




d(x, v) 





Making use of relations (30-32) one may reduce the terms on the right hand side of 
Eq.(33) to more simplified forms. Thus 



dP d dP 
V.F P(x,v,t) = -v— + P-(vP)-K(x)— , 



(34) 



a 



POO 

; V- / (FxV-r-F^x-^O) 
Jo 



d(x~ 



d(x, v) 



d 2 P d 2 P 
drP(x, v, t) = a 2 coTr^r + Q? c r 



dv 2 



dvdx 



(35) 



a 2 V- / r(F 1 V_ r -F 1 (^,0) 
Jo 



d(x T ,v T ) 



d(x, v) 



drP(x, v, t) 



-a 2 c 2 



d 2 P 

dv 2 



(36) 
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a 



PCX) 

! V- / r(F 1 V_ r -F 1 (^,OV^-Fo(x- T ,0) 
Jo 



d(x T , v T ) 



d(x, v) 



drP(x, v, t) 



a 2 Ci 



2 d 2 P + v d 3 P + ^ 9"? _ d 3 
<9t><9a; dv 2 dx dv 3 dv 3 



, (37) 



where 



£o = / °°(e(Oe(*-T)>dr 

5i = /o°°r(^)^-r)>dr [ ■ (38) 

£ 2 = /o°°r -r)>dr 
The final equation for the average motion corresponding to an adiabatic stochastic 

evolution in phase space is, 



dP „d , „, „, s dP ,, ,9 2 F 2 cfP 



+ a 2 5i 



d 3 P d 3 P d 3 

^^ + K(x) 1 --P—(vP) 



dv 2 dx 



dv 3 



dv 3 



(39) 



The remarkable departure from the standard form of Fokker-Planck equation is 
thus evident in Eq.(39) since it contains third derivative terms. The magnitude of 
their contribution is dependent on how much 'unperturbed' x and v vary during r 
which is of the order of \fi\~ 1 - We also point out that in the above derivation Bourret's 
decoupling approximation [10] has been used as in the treatment of linear equation in 
Sec.II. 
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V. Adiabatic fluctuations in velocity space 



We now consider the motion of a particle with velocity v in presence of fluctuations 
a£(t) which is adiabatically slow. The equation of motion is given by 



(40) 



The corresponding problem of a Brownian particle with fast fluctuations is a century- 
old problem in physical science, in general. Following the procedure described in the 
earlier section we first identify the perturbed and the unperturbed part of F, i.e., 



F = -[3v , Fi = 



(41) 



and calculate the Jacobian 



dv 



dv 



for the mapping v — > v T for the 'unperturbed' motion 



dv' 



dv 



/3r 



(42) 



Also note that 



V_ r = e 



-fir 



d_ 

dv 



(43) 



The evolution of the probability distribution function P(v, t) is then given by (terms 
of the order a 2 ) 



(44) 
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(45) 



where 

Ci 2 = Ci - C 2 

ci = / °°(^)e(*-T)>dr 

C2 = JS°T(t(t)t(t-T))dT 

While in the absence of the third term, the first two terms on the right hand side 
of Eq. (44) correspond to drift and diffusion terms in the Fokker-Planck description of 
an Ornstein-Uhlenbeck process, the third derivative term precludes the possibility of 
any straight-forward interpretation of the equation. Similar equations with third order 
noise, although not very common, however may be encountered [7] in the treatment 
of trimolecular reactions and also in quantum optics describing optical bistability in 
terms of associated Wigner distribution function for the reduced density operator in 
symmetrical ordering for the radiation field. 



We now return to the Eq.(44) which after some modification becomes 
3 n/ , n d r „, „ d 2 P(v,t) d 3 P(v,t) 



(46) 



a 2 c 3 



(47) 



where 

D± = a 2 (c 12 -3c 3 (3) 
D 2 = 

We now transform the Eq.(46) to Fourier space by defining the conditional proba- 
bility P(v,t\vo,0) and its Fourier transform as 

1 r +oc ; ~ 

(48) 



P(v, t\v , 0) = — / dke tkv P(k, t\v , 0) , 

ZTT J-oo 
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to obtain 



d ~ dP 

-P(k, t\v , 0) = -f3(k + D 2 k 3 )— - (D 1 + 3pD 2 )k 2 P 



(49) 



The linear partial differential equation (49) can be solved by the method of charac- 
teristics. For the initial condition (at time t — 0) 



P(v,0\v o ,0) =5(v-v ) , 



(50) 



the solution is 



P(k,t\v o ,0) 



(l + Bk 2 ) A 



exp 



-ikvn 



fit) 



1 + Bk 2 



where 



/(*) = 



-20t 



A = c 12 /2(3c 3 
and B = a 2 c 3 {l-f(t)} 
It is easy to check that Eq.(51) satisfies 



(51) 



(52) 



P*(k,t\v o ,0) = P(-k,t\v ,0) 



and the characteristic function (51) is of stable form. 

The conditional probability density P(v, t\vo, 0) is obtained by inverse Fourier trans- 
formation of Eq.(51) and is given by, 



1 r+ c 
P(v,t\v ,0) = —J 



oo dk (l + Bk^ eXP 



ikv — ikvo\ 



fit) 



1 + Bk 2 



(53) 
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It is evident that P(v,t\v ,0) is a Fourier transform of a stable characteristic func- 
tion. Hence the solution (53) forms a stable distribution in the variable v. The Eq.(53) 
is one of the important results of this paper. 

Although an explicit expression for P(v,t\v o ,0) is difficult to obtain, closed-form 
solutions for P(v,t\v o ,0) for stationary state can be easily obtained. In the long time 
limit the characteristic function (51) reduces to its asymptotic form 

which results a steady state distribution of stable form. Explicitly for small A, i.e., 
large (3 this is given by 

\ v \ A + l h 

Pst(v) = y e . (55) 

2 A D 2 2 F(A) v 2 

It is interesting to note that the dominant behavior of P st (v ) for small v. This is 
given by a power law of the form 



Pst(v) 



\v\' 



~ \v 



-l+A 



(56) 



Such power law behavior is also apparent for Levy processes [9,11] but for the large 
v regime. 

Although explicit solution for probability distribution P(v,t\v ,0) is difficult to ob- 
tain for arbitrary time, however, a few statistical properties of the process can be 
obtained from the calculation of variances and higher moments. For convenience, we 
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define such moments by subtracting the mean motion of the variables, i.e., we calculate 
the moments of At>(= v — t> e~ /3 *). Thus we write 



/+oo 
(v-v o e-^) m P{v,t\v o ,0) , 
-oo 



(57) 



or more explicitly 



1 r+°° 1 

(\Av\ m ) = — / dk- — - 

^ 1 1 27ri-oo (l + Bk 2 ) A 



exp 



After some algebra we get 



\Av\ m ) = H) 



+oo 

"' ' dk- 



l + Bk 2 ) A 



exp 



'ikv e-^^l-(l + Bk 2 ) ^} 

/+°° 
dAv\Av\ m e lkAv 
-oo 

ifcuoe-^jl- (l + Bk 2 )^} 



(58) 



d m 5{k) 
dk m 



(59) 



(60) 



In principle, using the property of Dirac 5-function and appropriate integrations any 
moment can be calculated from the above relation. We quote the results explicitly for 
the first three moments. 

For m=l (|Au|) = 
For m = 2 (\Av\ 2 ) = ^(1 - e - 2(St ) 
For m = 3 {\Av\ 3 } = 3a 2 c 3 v e^ l3t (l - e" 2/3t ) 
It is thus evident that unlike Levy processes [9,11] the moments are finite. 

We thus observe that because of the linear dissipation f3, a system driven by adiabatic 
noise reaches a steady state which is stable. However, since the noise is of external 
origin, the outward flow of energy due to linear dissipation cannot balance the inward 
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flow of energy supplied by the adiabatic fluctuations and hence a fluctuation-dissipation 
relation cannot be conceived in this case. 

VI. Conclusions 

In conclusion, we consider herein a dynamical system driven by external adiabatic 
fluctuations. Based on the 'adiabatic following approximation' we have made a sys- 
tematic separation of time-scales to carry out an expansion in a\fi\~ 1 to obtain a linear 
differential equation for the average solution, where a is the strength of fluctuation and 
is the largest eigenvalue of the unperturbed system. The main results of this study 
can be summarized as follows : 

(i) The probability distribution functions obey the differential equations of motion 
which contain third-order terms beyond the usual Fokker-Planck terms. The adiabatic 
fluctuations thus may give rise to non-Gaussian noise. 

(ii) We have examined in detail the corresponding equation in velocity space and 
the characteristic function is shown to obey a simple third-order differential equation 
which can be solved exactly in closed form. The characteristic function is found to be 
of stable form. 

(iii) Although third-order noise, in general, leads to serious interpretative difficulties 
in terms of truly probabilistic considerations in several cases, we show that in the 
present problem of adiabatic stochasticity in velocity space, statistical properties of the 
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processes are more transparent. It is of special interest to note that in contrast to Levy 
process all the variances and higher moments are finite and probability distribution is 
of stable form. 

(iv) Because of linear dissipation, the system driven by adiabatic fluctuations reaches 
a stable steady state. 

(v) For small arguments the probability distribution function obeys a power law 
behavior which is reminiscent of Levy processes. 

The stochastification by adding rapid fluctuating terms had been applied earlier to 
a wide variety of physical problems described by linear relaxation equations [15], hy- 
drodynamic equations [16], Maxwell equations in a medium [17], Boltzmann equation 
[18] etc. Our present analysis shows that the present method might reveal interest- 
ing consequences in such cases where the added fluctuating terms, in question, are 
adiabatically slow. We hope to address such issues in future communications. 
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